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f^*) , Abstract 

The problem of Subgraph Isomorphism is defined as follows: Given a pattern H and 
a host graph G on n vertices, does G contain a subgraph that is isomorphic to HI Eppstein 
[SODA 95, J'GAA 99] gives the first linear time algorithm for subgraph isomorphism for 
a fixed-size pattern, say of order k, and arbitrary planar host graph, improving upon the 
0(n v ^)-time algorithm when using the "Color-coding" technique of Alon et al [J'ACM 95]. 
Eppstein's algorithm runs in time kP^n, that is, the dependency on k is superexponential. 
We solve an open problem posed in Eppstein's paper and improve the running time to 2°^n, 
that is, single exponential in k while keeping the term in n linear. Next to deciding subgraph 
L—i • isomorphism, we can construct a solution and count all solutions in the same asymptotic 

running time. We may enumerate uj subgraphs with an additive term 0{ujk) in the running 
time of our algorithm. We introduce the technique of "embedded dynamic programming" 
on a suitably structured graph decomposition, which exploits the topology of the underlying 
embeddings of the subgraph pattern (rather than of the host graph). To achieve our results, 
we give an upper bound on the number of partial solutions in each dynamic programming 
step as a function of pattern size — as it turns out, for the planar subgraph isomorphism 
problem, that function is single exponential in the number of vertices in the pattern. 

^ ! 1 Introduction 

ON 

In the literature, we often find results on polynomial time or even linear time algorithms for 
NP-hard problems. Take for example the NP-complete problem of computing an optimal tree- 
decomposition of a graph. Bodlaender [3] gives an algorithm in time 0(n) for this problem — 
restricted to input graphs of constant treewidth. The Graph Minor Theory developed by Robert- 
son and Seymour implies amongst others that there is an 0(n 3 ) algorithm for the disjoint path 
problem, that is for finding disjoint paths between a constant number of terminals. Taking a 
closer look at such results, one notices that a function exponential in size of some constant c 
is hidden in the O-notation of the running time — here, c is the treewidth and the number of 
terminals, respectively. In another line of research, parameterized complexity, the primary goal 
is to rather find algorithms that minimize the exponential term of the running time. The first 
step here is to prove that such an algorithm with a separate exponential function exists, that 
is, that the studied problem is fixed parameter tractable (FPT) [131 HI EI] ■ Such problem has 
an algorithm with time complexity bounded by a function of the form f[k) ■ n°^ l \ where the 
parameter function f is a computable function only depending on k. The second step in the 
design of FPT-algorithms is to decrease the growth rate of the parameter function. 

We can identify two different trends in which running times of exact algorithms are improved. 
First, one can decrease the degree of the polynomial term in the asymptotic running time, and 
second, one can focus on obtaining parameter functions with better exponential growth. In 
the present work, we achieve both goals for the computational problem Planar Subgraph 
Isomorphism. 
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Subgraph Isomorphism generalizes many important graph problems, such as Hamil- 
tonicity, Longest Path, and Clique. It is known to be iVP-complete, even when restricted 
to planar graphs [18J. Until now, the best known algorithm to solve SUBGRAPH ISOMORPHISM, 
that is to find a subgraph of a given host graph isomorphic to a pattern H of order k (the 
number of vertices in H), is the naive exhaustive search algorithm with running time 0(n k ) 
and no FPT-algorithm can be expected here [13]. For a pattern H of treewidth at most t, Alon 
et al. p] give an algorithm of running time 2°^n°^ . For PLANAR SUBGRAPH ISOMORPHISM, 
given planar pattern and input graph, some considerable improvements have been made mostly 
during the 90's. The first improvement was provided by Plehn and Voigt [22], with running 
time 0(k k )n (^\ Using the elegant Color-coding technique of Alon et al. p], one can devise an 
algorithm of running time 2 ^n°( y ^ k ^ > . The current benchmark has been set by Eppstein [14] to 
k°^n, by employing graph decomposition methods, similar to the Baker-approach [2] for ap- 
proximating NP-complete problems on planar graphs. Eppstein's algorithm is actually the first 
FPT-algorithm for Planar Subgraph Isomorphism with k as parameter. Eppstein poses 
three open problems: a) whether one can extend the technique in [l] to improve the dependence 
on the size of the pattern from k°^ to 2°^ for the decision problem of subgraph isomorphism; 
and whether one can achieve similar improvements b) for the counting version and c) for the 
listing version of the subgraph isomorphism problem. 

Our results. In this work, we do not only achieve this single exponential behavior in k for all 
three problems — without applying the randomized coloring technique — we also keep the term in 
n linear. That is, we give an algorithm for Planar Subgraph Isomorphism for a pattern H of 
order k with running time 2°^n. Next to deciding subgraph isomorphism, we can construct a 
solution and count all solutions in the same asymptotic running time. We may list u subgraphs 
with an additive term 0(uik) in the running time of our algorithm. Our algorithm also improves 
the time complexity of the previous approach [17] for patterns of size k £ o(y / nlogn). 

The novelty of our result comes from embedded dynamic programming, a technique we find 
interesting on its own. Here, one decomposes the graph by separating it into induced subgraphs. 
In the dynamic programming step, one computes partial solutions for the separated subgraphs, 
that are updated to an overall solution for the whole graph. In ordinary dynamic programming, 
one would argue how the subgraph pattern hits separators of the host graph. Instead, in 
embedded dynamic programming for subgraph isomorphism, we proceed exactly the opposite 
way: we look at how separators can be routed through the subgraph pattern. As a consequence, 
we bound the number of partial solutions not by a function of the separator size of the host graph, 
but by a function of the pattern size — as it turns out, for the planar subgraph isomorphism 
problem, that function is single exponential in the number of vertices of the pattern. To obtain a 
good bound on the parameter function, we apply several fundamental enumerative combinatorics 
results in the technical sections of this work. Next to the number of cycles and face-vertex 
sequences in embedded graphs, these counting results give upper bounds on the number of 
planar triangulations and planar embeddings of the pattern. 

Our algorithm is divided into two parts with the second part being the aforementioned 
embedded dynamic programming. For keeping the time complexity of our algorithm linear in the 
size of the host graph, we give a fast method for computing a graph decomposition with special 
properties: Sphere-cut decompositions are natural extensions of tree-decompositions to plane 
graphs, where the separator vertices are connected by a Jordan curve. In embedded dynamic 
programming we use sphere-cut decompositions with separators of size linearly bounded by the 
size of the subgraph pattern. 

Theorem 1.1 Let G be a planar graph on n vertices and H a pattern of order k. We can 
decide if there is a subgraph of G that is isomorphic to H in time 2°^n. We find subgraphs 
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and count subgraphs of G isomorphic to H in time 2°^n and enumerate uj subgraphs in time 

It is worth mentioning that for /c-Longest Path on planar graphs, the authors of [12] 
give the first algorithm with time complexity sub exponential in the parameter value. The 
algorithm has running time 2° lyy ^ l n + 0(n 3 ), employing the techniques Bidimensionality and 
topology-exploiting dynamic programming. Bidimensionality Theory employs results of Graph 
Minor Theory by Robertson and Seymour for planar graphs |23j and other structural graph 
classes to algorithmic graph theory (entry [6] , for a survey [7] ) . Unfortunately, Bidimensionality 
does only work for finding specific patterns in a graph, such as /c-paths, but not for subgraph 
isomorphism problems in general. For a survey on other planar subgraph isomorphism problems 
with restricted patterns, please consider [14] . 

Organization. After giving some definitions in Section [21 we show in Section [3] how to ob- 
tain a sphere-cut decomposition of small width. In Section [5] we restrict Planar Subgraph 
Isomorphism to Plane Subgraph Isomorphism. We first give some technical lemmas in 
Section 14.11 to bound the number of ways a separator of the sphere-cut decomposition can be 
routed through a plane pattern. We describe and analyze embedded dynamic programming in 
Section 14.21 followed by subsuming the entire algorithm for Plane Subgraph Isomorphism 
in Section [4.31 In Section [5] we bound the number of drawings of the pattern and show how to 
solve Planar Subgraph Isomorphism. 

2 Preliminaries 

Subgraph isomorphism. Let G, H be two graphs. We call G and H isomorphic if there 
exists a bijection v : V{G) -> V(H) with {v,w} G E(G) ^ {v(v),v(w)} G E(H). We call H 
subgraph isomorphic to G if there is a subgraph H' of G isomorphic to H. 

Branch Decompositions. A branch decomposition (T, fi) of a graph G consists of an un- 
rooted ternary tree T (i.e., all internal vertices have degree three) and a bijection \i : L — > E(G) 
from the set L of leaves of T to the edge set of G. We define for every edge e of T the middle 
set mid(e) C V(G) as follows: Let T\ and T2 be the two connected components of T \ {e}. 
Then let Gi be the graph induced by the edge set {//(/) : f £ Ln V(Tj)} for i G {1, 2}. The 
middle set is the intersection of the vertex sets of G\ and G2, i.e., mid(e) := V(G\) n ViG?). 
The width bw of (T,fj) is the maximum order of the middle sets over all edges of T, i.e., 
bw((T, fj,)) := max{| mid(e)| : e G T}. An optimal branch decomposition of G is defined by a 
tree T and a bijection [i which together provide the minimum width, the branchwidth bw(G). 

Plane graphs and equivalent embeddings. Let £ be the unit sphere. A plane drawing 
or planar embedding of a graph G with vertex set V(G) and edge set E(G) maps vertices to 
points in the sphere, and edges to simple curves between their end vertices, such that edges do 
not cross, except in common end vertices. A plane graph is a graph G together with a plane 
drawing. A planar graph is a graph that admits a plane drawing. For details, see e.g. [10] . 
The set of faces F(G) of a plane graph G is defined as the union of the connected regions of 
£ \ G. A subgraph of a plane graph G, induced by the vertices and edges incident to a face 
/ G F(G), is called a bound of /. If G is 2-connected, each bound of a face is a cycle. We 
call this cycle face-cycle (for further reading, see e.g. [10J ) . For a subgraph H of a plane graph 
G, we refer to the drawing of G reduced to the vertices and edges of H as a subdrawing of G. 
Consider any two drawings G\ and G2 of a planar graph G. A homeomorphism of G\ onto G2 
is a homeomorphism of E onto itself which maps vertices, edges, and faces of G\ onto vertices, 
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edges, and faces of G2, respectively. We call two planar drawings of the same graph equivalent, 
if they are homeomorphic. 

Theorem 2.1 (e.g. fH$) Every 3-connected planar graph has a unique embedding in a sphere 
£ up to homeomorphism. 

Triangulations. We call a plane graph G a planar triangulation or simply a triangulation if 
every face in F(G) is bounded by a triangle (a cycle of length three). If H is a subdrawing of a 
triangulation G, we call G a triangulation of H. 

Nooses and combinatorial nooses. A noose of a S-plane graph G is a simple closed curve 
in £ that meets G only in vertices. Prom the Jordan Curve Theorem, it then follows that 
nooses separate S into two regions. Let V(N) = N n V(G) be the vertices and F(N) be the 
faces intersected by a noose N. The length of N is the number |V(iV)| of vertices in V(N). 
The clockwise order in which N meets the vertices of V(N) is a cyclic permutation tt on the set 
V(N). 

Remark 2.2 Let a plane graph H be a subdrawing of a plane graph G. Every noose N in G is 
also a noose in H and Vh{N) C Vg(N). 

A combinatorial noose Nc = [vq, fo, v\, f\, . . . , fk-i, Ufe] in a plane graph G is an alternating 
sequence of vertices and faces of G, such that 

• fi is a face incident to both Vi, fj+i for all i < k, 

• t> = Vk and the vertices v\ , . . . , Vk are mutually distinct and 

• if fi = fj for any i ^ j and i,j = 0,...,k—l, then the vertices Vi, i>i+i, Uj, and Uj+i do not 
appear in the order (vi,Vj,Vi+i,Vj+i) on the bound of face /, = fj. 

The length of a combinatorial noose [vo, fo, v\, f\, . . . , fk-i,Vk] is k. 

Remark 2.3 The order in which a noose N intersects the faces F(N) and the vertices V(N) 
of a plane graph G gives a unique alternating face-vertex sequence of F(N) U V(N) which is a 
combinatorial noose Nq- Conversely, for every combinatorial noose Nc there exists a noose N 
with face-vertex sequence Nc- 

We may view combinatorial nooses as equivalence classes of nooses, that can be represented by 
the same face-vertex-sequence. 

Sphere cut decompositions. For a S-plane graph G, we define a sphere cut decomposition 
or sc- decomposition (T, fi, tt) as a branch decomposition which for every edge e of T has a 
noose iV e that cuts E into two regions Ai and A2 such that Gi C Aj U -/V e , where Gi is the 
graph induced by the edge set {fi(f) : / G L n V(Ti)} for i 6 {1, 2} and TiUT 2 = T\ {e}. Thus 
A^e meets G only in y(iV e ) = mid(e) and its length is | mid(e)|. The vertices of every middle 
set mid(e) = V{G\) n V{G2) are enumerated according to a cyclic permutation tt on mid(e). 

The following two propositions will be crucial in that they give us upper bounds on the 
number of partial solutions we will compute in our dynamic programming approach. With both 
propositions, we will bound the number of combinatorial nooses in a plane graph by the number 
of cycles in the triangulation of some auxiliary graph. With the second proposition we bound 
the number of non-equivalent embeddings of planar graphs. 

Proposition 2.4 (J^j) No planar n-vertex graph has more than 2 L53n simple cycles. 

Proposition 2.5 (]26$) The number of non-isomorphic maximal planar graphs on n vertices 
is approximately 2 3 - 24n . 
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Proposition 12.51 also gives a bound on the number of non-isomorphic triangulations. Any em- 
bedding of a maximal planar graph G must be a triangulation, otherwise G would not be 
maximal. With Theorem 12. 1\ every maximal planar graph has a unique embedding which is a 
triangulation. On the other hand, every triangulated graph is maximal planar. 

3 Computing sphere-cut decompositions in linear time 

In this section we introduce an algorithm for computing sc-decompositions of bounded width. 
Let H be a connected subgraph of G with |V(i?)| = k, and let v £ V(H). Then H is a 
subgraph of the induced subgraph G v of G, where G v = G[S] with S = {w £ S | dist(u, w) < k] 
(dist(f , w ) denotes the length of a shortest path between v and w m G). This observation helps 
us to shrink the search space of our algorithm by cutting out chunks of G of bounded width 
and solve subgraph isomorphism separately on each chunk. With the algorithm of Tamaki |25j . 
one can compute a branch decomposition of G v of width < 2k + 1, following similar ideas as in 
the approach of Baker [2] for tree decompositions. With some simple modifications, we achieve 
the same result for sc-decompositions. In Appendix [A] we prove the following lemma and give 
an algorithm that computes a sc-decomposition of bounded width in linear time. 

Lemma 3.1 Let G be a plane graph with a rooted spanning tree whose root-leaf-paths have 
length at most k. We can find an sc-decomposition of width 2k + 1 in time 0{kn). 

4 Plane subgraph isomorphism 

In this section, we study the subgraph isomorphism problem on patterns and host graphs that 
are embedded in a sphere S. In Section [5] we carry over our results to planar graphs. We 
first introduce some topological tools that allow us to define a refined dynamic programming 
approach. At every step of the dynamic programming approach, we compute all possibilities of 
how a combinatorial noose N corresponding to a middle set of the sc-decomposition (T, fi, tt) 
of G can intersect a subdrawing equivalent to pattern H. Each intersection gives rise to a 
combinatorial noose of H. See Figure Q] for an illustration. 



Figure 1: On the left, we have a plane graph G with an emphasized subdrawing H intersected by 
a combinatorial noose N indicated by dashed lines. On the right, we have the same graph G with a 
different copy of H intersected by N. 

The running time of the algorithm crucially depends on the number of combinatorial nooses 
in H . The aim of this section is to prove the following: 

Theorem 4.1 Let G be a plane graph on n vertices and H be a plane graph on k <n vertices. 
We can decide if there is a subdrawing of G that is equivalent to H in time 2°^n. We can find 
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and count subdrawings equivalent to H in time 2°^n, and enumerate uj subdrawings in time 
2°^n + 0{ujk). 

4.1 Combinatorial nooses in plane graphs 

For a refined algorithm analysis we now take a close look at combinatorial nooses of plane 
graphs. In particular we are interested in counting the number of combinatorial nooses. In this 
subsection, we will prove the following lemma: 

Lemma 4.2 Every plane k-vertex graph has 2°^ combinatorial nooses. 

Before proving this lemma, we show that every combinatorial noose of a plane graph on k 
vertices corresponds to a cycle in some other plane graph on at most 0(k) vertices. First we 
relate combinatorial nooses in a planar triangulation H' to the cycles in H' . In a second step 
we relate combinatorial nooses of a 3-connected plane graph H to cycles in the triangulations of 
H. Finally, we will show that for any plane graph H there is an auxiliary graph H*, such that 
the combinatorial nooses of H can be injectively mapped to the cycles of the triangulations of 
H* . From Proposition 12.41 we know an upper bound on the number of cycles in planar graphs, 
which we employ to prove Lemma l4.2i 

Lemma 4.3 Let H be a planar triangulation and Nc = [vq, fo, V\, f\, . . . , fk-\, v k\ a combina- 
torial noose of H. Then for every pair of consecutive vertices Vi,Vi+\ in Nc, there is an edge 
Vi,Vi + \ in E(H). That is, the sequence [vo,vi, . . . , v^] is a simple cycle in H if\V(Nc)\ > 2, 
and if \V(Nc)\ = 2, it corresponds to a single edge in H. 

Proof. Since H is triangulated, we have that every fi £ Nq is bounded by a triangle A where 
Vi,Vi+\ are two of the three vertices of A and Vi, Uj+i have an edge in common. Since vertices 
occur only once in Nc, fi is unique in Nc if |V(iVc)| > 2, that is, there is no fj £ Nc with 
i / j and fi = fj. Hence we map each fi one-to-one to edge &i = {vi,Vi + ±} and get a cycle 
[vo, eo, vi, e±, . . . , e^-i, For ^(./VcOl = 2, /o and f\ are incident faces to edge {vq,v{\. For 
an illustration, see Figure □ 




Figure 2: On the left, there is a triangulation with a combinatorial noose Nc indicated by dashed lines. 
On the right, we have mapped the noose to a cycle indicated by dashed lines. 

Lemma 4.4 Let H be a 3-connected plane graph and Nc = [vo, fo, Vi, fx, . . . , fk-i, Vk\ a combi- 
natorial noose of H with \V(Nq)\ > 2. Then there exists a planar triangulation H' of H, such 
that [vo,vi, . . . , Vk] is a cycle in H' . 

Proof. We proceed in two phases. First we iteratively add edges to H and transform Nc into 
another combinatorial noose such that every two consecutive vertices in Nc have a common 
edge. Then we triangulate the resulting graph. 

For every pair of consecutive vertices Vi, Uj+i in Nc, if Vi,Vi+\ have no edge in common, 
add &i = {vi,Vi + \} to E{H). Thereby the drawing of splits fi into two new faces ff and 
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ff, bounded by face-cycle C a and C b respectively, where C a n C b = ej. Since Nc corresponds 
to a noose by Remark 12.31 and nooses are not self-intersecting, we observe the following for 
| V(iVc?) | > 2: for every fj = /j in Nc with j / i we have that both Vj,Vj + ± are in one of C a 
and C&. Thus, adding edge ej will not cross any other edge added in this process. In Nc, we 
replace fi by one of /? and and every /j = /j by /? if is bounded by C a and by f\ 
otherwise. Once we have an edge for every pair of consecutive vertices in Nc, we note that 
for every sub-sequence [vi, fi, Uj+i] of Nc the edge ej = is incident to face fi since, 

by 3-connectivity, edge is uniquely embedded in H. We then add edges arbitrarily to obtain 
a triangulation H' . By Lemma 14.31 the vertices of Nc correspond to a cycle in H' . For an 
illustration, see Figure El □ 



Figure 3: One the left, we see a face of our plane graph with a combinatorial noose Nc indicated with 
dashed lines. To the right, we have mapped the noose to a cycle indicated with dashed lines and dotted 
lines indicating the face triangulation. 

If H is not 3-connected, a problem may occur in the last step of the previous proof when 
triangulating H. Consider a sub-sequence [uj, fi, Uj+i] in Nc. We assume there already exists 
an edge ej = {vi,Vi + i} and Vi,Vi + \ separate H, that is, H is 2-connected. Then it may be the 
case that e« is not incident to fi, and thus, any triangulation of H has an edge crossing Nc- We 
surpass this problem in the general case by triangulating some auxiliary graph instead. For an 
edge e = {v, w} of a graph H we subdivide e by adding a vertex u to V(H) and replacing e by 
two new edges e\ = {v,u} and e<i = {u, w}. In a drawing of H, we place point u in the middle 
of the drawing of e partitioning e into e\ and e^. 

Lemma 4.5 Let H be plane graph and Nc = [vq, fo, V\, fx, . . . , fk~\, v k] a combinatorial noose 
of H with \V(Nc)\ > 2. Let H* be obtained by subdividing every edge in E(H). There exists a 
planar triangulation H' of H* such that [vq,v\, ... ,Vk] is a cycle in H' . 

Proof. The combinatorial noose Nc is a combinatorial noose in H* , too. As for any two 
consecutive vertices t^Vj+i in Nc there is no edge in H* and each vertex in Nc is unique, we 
may add edges to H* as in the proof of Lemma 14.41 and triangulate H* . □ 

Proof of Lemma l4.2l If H is triangulated, we have with Lemma [4.3l that every combinatorial 
noose corresponds to a unique cycle in H. By Proposition 12.41 the- number of cycles in H is 
bounded by 2 1,53fc . Since for every edge of a cycle in H, we have two choices for a combinatorial 
noose to visit an incident face, we get the overall upper bound of 2 2,53fc on the number of 
combinatorial nooses. If H is plane, we have to count the triangulations either of H (Lemma l4.4p 
or of H* (Lemma I4.5p . By Proposition 12.51 and the comments below it, there are at most 2 3 24 ^ 
non- isomorphic triangulations on £ vertices. Let us denote this set of triangulated graphs by 
<!>. We note that H (resp. H*) is a subgraph of some graph of <£, say of all graphs in C 
with > !• Since every triangulated graph is 3-connected, we have with Theorem 12.11 that 
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every graph H' in &h has a unique embedding in E up to homeomorphism. The plane graph 
H (resp. H*) is then a subdrawing of a drawing equivalent to an arbitrary plane embedding 
of H' in S. Thus, the number of triangulations times the number of combinatorial nooses in 
each triangulation is an upper bound on the number of combinatorial nooses in H, here 2 5,77fc 
(resp. in H* , here 2 9 ' m ). □ 

For embedded dynamic programming on a sc-decomposition (T,fj,,w), we can argue with 
Remark 12.21 that if H is a subdrawing of G, then noose iV formed by the middle set mid(e) is 
a noose of H, too. Recalling Remark 12.31 the alternating sequence of vertices and faces of H 
visited by N forms a combinatorial noose Nq in H. 

This observation allows us to discuss the results from a combinatorial point of view without 
the underlying topological arguments. Instead of nooses we will refer to combinatorial nooses 
in the remaining section. 



4.2 Embedded dynamic programming 

In embedded dynamic programming, the basic difference to usual dynamic programming is 
that we do not check for every partial solution for a given problem if or how it lies in the graph 
processed so far. Instead, we check how the graph that we have processed so far is intersecting the 
entire solution, that is how the graph is embedded into our solution. For subgraph isomorphism, 
we compute every possible way the processed subdrawing G su b of G is embedded in the plane 
pattern H up to homeomorphism, subject to how the bound of G su b intersects H. This bound 
is a combinatorial noose N separating G su b from the rest of G. The number of solutions we get 
is bounded by the number of combinatorial nooses in H we can map N onto. We describe the 
algorithm in what follows. 

Dynamic programming. We root sc-decomposition (T,fi,ir) at some node r G V(T). For 
each edge e € T, let L e be the set of leaves of the subtree rooted at e. The subgraph G e of G 
is induced by the edge set {fJ>(v) \ v £ L e }. The vertices of mid(e) form a combinatorial noose 
N that separates G e from the residual graph. 

Assuming H is a subdrawing of G, the basic idea of embedded dynamic programming is 
that we are interested in how the vertices of the combinatorial noose N are intersecting faces 
and vertices of H. Since every noose in G is a noose in H, we can map N to a combinatorial 
noose N H of H, bounding (clockwise) a unique subgraph H su b of H. 

In each step of the algorithm, all solutions for a sub-problem in G e are computed, namely 
all possibilities of how N is mapped onto a combinatorial noose N H in H that separates H su b 
from the rest of H, where H su b C H is isomorphic to subgraphs of G e . For every middle set, 
we store this information in an array. It is updated in a bottom-up process starting at the 
leaves of (T, /i,7r). During this updating process it is guaranteed that the 'local' solutions for 
each subgraph associated with a middle set of the sc-decomposition are combined into a 'global' 
solution for the overall graph G. 

Step 0: Initializing the middle sets. Let G be a plane graph with a rooted sc- 
decomposition (T, fi, 7r) and let H be a plane pattern. For every middle set mid(e) of (T, (a, tt) 
let iV be the associated combinatorial noose in G with face- vertex sequence of F(N) U V(N). 
Let £ denote the set of all combinatorial nooses of H whose length is at most the length of 
N. We now want to map N order preserving to each N H E £. We map vertices of N to 
both vertices and faces of H. Therefore, we consider partitions of V(N) = Vi(N)0V2(N) where 
vertices in V\(N) are mapped to vertices of V(H) and vertices in V^(A^) to faces of F{H). We 
define a mapping 7 : V(N) U F(N) — > V{H) U F{H) relating iV to the combinatorial nooses 
in £. For every N H G £ on faces and vertices of set F(N H ) U V(N H ) and for every partition 
Vi(N)0V 2 (N) of V(N) mapping 7 is valid if 
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a) 7 restricted to V\(N) is a bijection to V(N H ); 

b) for every v G Vi(iV) we have 7(f) G V(N H ), and for every v G ^(-/V) we have 7(1?) G F(N H ); 

c) for every / G F(iV) we have 7(/) G F(N H ); 

d) for every pair UftjVj G V^(iV) such that [7(11^), /, j(vj)] is a subsequence of iV^ for a face 
/ G F(N H ) and for every vertex m G V(iV) with lying inbetween Vh and in the sequence 
N, we have j(v{) = /; 

e) for every v-i G V(JV) and subsequence [fi-i,Vi, fi] of iV, if j(vi) G F(N H ), we have 7(/j-i) = 
tK) = 7(/i); 

/) for every pair of vertices Wi,Wj in V(iV^): if {w^iOj} G E{H) then {7 -1 (itfi), 7 _1 ( w j)} G 
S(G). 

Items a) to c) say where to map the faces and vertices of N to. Items d) and e) make sure that 
if two vertices Vh,Vj in sequence N = [. . . , Vh, ■ ■ -,Vj, ■ ■ ■] are mapped to two vertices Wi,Wi + \ 
that appear in sequence N H as [. . . , w%, fi,Wi+i, ■ ■ .] then every face and vertex inbetween Vh, Vj 
in sequence N (here underlined) is mapped to face fi. Item /) rules out the invalid solutions, 
that is, we do not map a pair of vertices in G that have no edge in common to the endpoints of 
an edge in H. We do so because if H is a subdrawing of G then an edge in H is an edge in G, 
too. For an illustration, see Figure |H 




Figure 4: On the left, we have a plane graph G with a subdrawing H emphasized. A combinatorial 
noose N separating subgraph G e is indicated by dashed lines. The vertices of N are full and empty 
circles and the faces triangles. In the middle, we have H and indicate to which faces (big triangles) of H 
vertices and faces of N are mapped by 7. This gives us combinatorial noose N H on the right, separating 
subgraph H sub . 

We assign an array A e to each mid(e) consisting of all tuples (N H ,j(N)} each representing 
a valid mapping j(N) from combinatorial noose iV corresponding to mid(e) to a combinatorial 
noose N H G £. The vertices and faces of N are oriented clockwise around G e . Without loss of 
generality, we assume for every (N H ,j(N)) G A e the orientation of N H to be clockwise around 
the subgraph H su b of H isomorphic to a subgraph of G e . 

Step 1: Update process. We update the arrays of the middle sets in post-order manner 
from the leaves of T to root r. In each dynamic programming step, we compare the arrays of 
two middle sets mid(e), mid(/) in order to create a new array assigned to the middle set mid(g), 
where e, / and g have a vertex of T in common. From [12] we know about a special property 
of sc-decompositions: namely that the combinatorial noose N g is formed by the symmetric 
difference of the combinatorial nooses N e ,Nf and that G g = G e U Gf. In other words, we are 
ensured that if two solutions on G e and Gf bounded by N e and Nf fit together, then they form 
a new solution on G g bounded by N g . We now determine when two solutions represented as 
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tuples in the arrays A e and Af fit together. We update two tuples (N^ ,j e (N e )) G A e and 
(Nf ,7/(iV/)) G Af to a new tuple in A g if 

• for all v G V(JV e ) n V(N f ), 7e (u) = 7/ (<w); 
. for all / G F(iVe) D F(JV/), 7e (/) = 7/ (/); 

• for the subgraph of -ff separated by A^ and the subgraph ffy of separated by A^, we 
have that E(H e ) n = and V(£T e ) D C { 7 (u) | u G V(N e ) n V(iV»}. 

If N e and Nf fit together, we get a valid mapping 7s : iV a — » A^ as follows: 

• for every x G ( V(JV e ) U F{N e ) ) n ( V(Nf) U F(Nf) ) n ( V(JV fl ) U F(N g ),) we have 7e (x) = 

7/0*0 = 7ff(«); 

. for every y G ( V(N e ) U F(N e ) ) \ ( V(Nf) U F(JV» ) we have 7e (y) = lg (y); 

• for every 2 G ( V(N f ) U F(Nf) ) \ ( F(AT e ) U F(N e ) ) we have 7/ (z) = j g (z). 

We have that 7s is a valid mapping from N g to the combinatorial noose N g that bounds 
subgraph H g = H e U iJy. Thus, we add tuple (Nf ,j g (N g )) to array 

Step 2: End of DP If, at some step, we have a solution where the entire subgraph H is 
formed, we exit the algorithm confirming. That is, if H = H e U Hf and Hi is bounded by ATj 
(for both i G {e, /}) then the combinatorial noose A^ is bounding the subgraph of G isomorphic 
to H. We are able to output this subgraph by reconstructing the solution top-down in (T, fj,, it). 
If at root r no subgraph isomorphic to H has been found, we output 'FALSE'. 

Correctness of DP Let plane graph H be a subdrawing of G. We have seen already in Step 
how we map every combinatorial noose of G that identifies a separation of G via a valid mapping 
7 to a combinatorial noose of H determining a separation of H. Every edge of H is bounded by 
a combinatorial noose A^ of length two, which is determined by tuple (A r ^, 7 (A r )} in an array 
assigned to a leaf edge of T. We need to show that Step 1 computes a valid solution for N g from 
A^e and Nf for incident edges e, /, g. We note that the property that the symmetric difference 
of the combinatorial nooses N e and Nf forms a new combinatorial noose A^ is passed on to 
the combinatorial nooses N^,Njr and Njf of H, too. If the two solutions fit together, then 
H e of H separated by A^ and subgraph Hf of H separated by A^ only intersect in the image 
of V^Ae) D V(Nf). We may observe that N^ 1 and A 7 ^ intersect in a continuous alternating 
subsequence with order reversed to each other, i.e., N^ |Ar e nAf/ = N? \N e nN f , where A^ means 
the reversed sequence A^. Since every oriented A^ identifies uniquely a separation of E(H), 
we can easily determine if two tuples (N^ , 7e (JV e )) G A e and (Nj 1 ,jf(Nf)} G Af fit together 
and form a new subgraph of H. If H is a subdrawing of G, then at some step we will enter Step 
2 and produce the entire H. 

Running time analysis. We first give an upper bound on the size of each array. The number 
of combinatorial nooses in £ we are considering is bounded by the total number of combinatorial 
nooses in H, which is 2°^'^'^ by Lemma 14.21 The number of partitions of vertices of any 
combinatorial noose A^ is bounded by 2'^^' . Since the order of both A^ and A^ is given we only 
have 2|V(jy)| possibilities to map vertices of A^ to A^, once the vertices of Af are partitioned. 
Thus, in an array A e we may have up to 2°^ V(H ^ -2^^ ■ \V{H)\ tuples (N^,j(N e )). We first 
create all tuples in the arrays assigned to the leaves. Since middle sets of leaves only consist 
of an edge in G, we get arrays of size 0(\V(H)\ 2 ) which we compute in the same asymptotic 
running time. When updating middle sets mid(e), mid(/), we compare every tuple of one array 
A e to every tuple in array Af to check if two tuples fit together. We can compute the unique 
subgraph H e (resp. Hf) described by a tuple in A e (resp. Af), compare two tuples in A e ,Af 
and create a new tuple in A g in time linear in the order of V^(A^) and V(H). Since the size of 
A g is bounded by 2°^ H ^ ■ 2°^ N ^ , the update process for two middle sets takes the same 



10 



asymptotic time. Assuming sc-decomposition (T,fi,ir) of G has width u and |V(i7)| < u>, we 
get the following result. 

Lemma 4.6 For a plane graph G with a given sc-decomposition {T, fi, it) of G of width w and 
a plane pattern H on k < w vertices we can search for a subdrawing of G equivalent to H in 
time 2°M • n. 



4.3 The algorithm 

We present the overall algorithm for solving Plane Subgraph Isomorphism with running 
time stated in Theorem 14.11 

Algorithm 4.1: Plane subgraph isomorphism: PLSI. 
Input : Plane graph G; Plane pattern H of order k. 

1 Choose an arbitrary vertex v in G. 

2 Partition V(G) into So U Si U . . . U S e with Si = {w G V(G) : d\st(v, w) = i} 

3 for every Gi = G[Si U . . . U S i+ k] with 0<i<£-kdo 

4 Compute sc-decomposition (T, /U,7r) of Gi. 

5 Do embedded dynamic programming on (T, fi,n) to find a subgraph of Gi isomorphic 
to H and intersecting Si. 



Partitioning the vertex set in Line [TJ of Algorithm 14.11 PLSI, is a similar approach to the 
well-known Baker-approach [2]. Every vertex set S% contains the vertices of distance i to the 
chosen vertex v. Sq = {v} and £ is the maximum distance in G from v. The graph Gi in Line [I] 
is induced by the sets Si, ... , Si + k- As in [14], we may argue that every vertex in G appears in at 
most k subgraphs Gi. This keeps our running time linear in n. We can apply Lemma f3.1l to each 
Gi in Line[T]to a compute sc-decomposition (T, fj,, n) of width < 2k, by adding a root vertex r for 
the BFS tree and make r adjacent to every vertex in Sj. The dynamic programming approach 
can easily be turned into an algorithm counting subgraph isomorphisms (similar to |14j), by 
using a counter in the dynamic programming. Using an inductive argument, for every subgraph 
Gi in Line [1] we only compute subgraphs intersecting with vertices in S« and thus omit double- 
counting. We can also adopt our technique to list the subgraphs of G isomorphic to H. 



5 Planar subgraph isomorphism 

Now we consider the case when both pattern H and host graph G are planar but not embedded. 
However, we observe that if H is isomorphic to a subgraph of G, then for every planar embedding 
of G there exists a drawing of H that is equivalent to a subdrawing of G. Hence, we may 
simply embed G planarly, and run the algorithm of the previous section for all non-equivalent 
embeddings of H. The following lemma tells us that the number of times we call the algorithm 
is restricted, too. 

Lemma 5.1 Every planar k-vertex graph has 2°^ non- equivalent embeddings in S. 

Proof. By Proposition 12.5^ there are at most 2 3,24fc non-isomorphic maximal planar graphs on 
k vertices. Every planar graph H is a subgraph of a maximal planar graph. Every maximal 
planar graph has a unique embedding which is a triangulation. Thus, every embedding of H is 
a subdrawing of a triangulation of H. The number of such subgraphs is bounded by the number 
of edge subsets of H' , since for every edge subset of S C E(H') of same cardinality as E(H), 
H'[S] may be isomorphic to H. In this case, H'[S] then gives a possible embedding of H in E. 
Hence, the number of embeddings of H in S up to homeomorphism is bounded by 2 6,24fc . □ 
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Algorithm 5.1: Planar subgraph isomorphism. 



Input : Planar graph G, Planar pattern H of size k. 
Compute a planar embedding of G. 

if H triangulated or 3-connected then Return PLSI(G, H). 
for every non- equivalent embedding I of H do 
L Return PLSI(G,/). 



The whole algorithm We compute in Algorithm 15.11 every non-equivalent embedding of H 
using the constructive proof of Lemma 15.11 That is, we compute the set Ti of non-isomorphic 
maximal planar graphs in time proportional to its size using the algorithm in [20]. For every 
graph H' E Ti and every subdrawing / of H' we check whether / is isomorphic to H by using 
the linear time algorithm for planar graph isomorphism in [19^1 . By Lemma l5.lt we then call 
Algorithm 14.11 2°^ times, for each plane drawing / isomorphic to H. This ensures us that 
Algorithm 15.11 has running time as stated in Theorem 11.11 

6 Conclusion 

We have shown how to use topological graph theory to improve the results on the already 
mentioned variations of Planar Subgraph Isomorphism, solving the open problems posed 
in [14] and [12] . With the results of [15] , [14] extends the feasible graph class from planar graphs 
to apex-minor-free graphs. This cannot be done with the tools presented here. However, |11] 
devise a truly sub exponential algorithm for /c-Longest Path in //-minor-free graphs and thus 
apex-minor-free graphs, employing the structural theorem of Robertson and Seymour [2l] and 
the results of [HJ[5l[9]. Can the structure of H- minor- free graphs, be exploited for our purposes? 

It seems unlikely that our work can be extended to obtain a subexponential algorithm. 
The first reason, mentioned in the introduction, is that Bidimensionality applies to subgraphs 
with minor properties rather than to general subgraphs. Secondly, our enumerative bounds are 
either tight or of lower bound 2^( fc ) . We want to pose the open problem: Is Plane Subgraph 
Isomorphism solvable in time 2 o(7c )n°( 1 )? 

Acknowledgments. The author thanks Paul Bonsma, Holger Dell and Fedor Fomin for dis- 
cussions and comments of great value to the presentation of these results. 
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A SC-Decompositions in linear time 

For a plane graph G we define a radial graph Rg as follows: Rg is a bipartite graph with the 
bipartition F(G) U V(G). A vertex v G V{G) is adjacent in Rg to a vertex / G F(G) if and 
only if the vertex v is incident to the face / in the drawing of G. 

Let G be a plane graph with some vertex r G V{G) and Rg its radial graph. Let T be a 
spanning tree of G rooted at / that is determined by breadth first search. Choose a face / 
adjacent to root r. If the longest path from r to a leaf of T is I then the distance df in the 
radial graph Rq from vertex / to any other (face) vertex x is at most 2£ + 1. This is due to 
the fact that there exists an edge {/, r} in Rq, and for every edge in T there is a detour in Rg 
of at most two edges. [25] show how to obtain a branch decomposition of width df out of a 
BFS spanning tree rooted at r of the radial graplJl. Set / to be the outer face of G. Let T be 
a BFS spanning tree of Rq rooted at / and let I be the maximum distance in T from r to a 
leaf. We give now a compact presentation of the algorithm of [25] and show that it translates 
to constructing a sphere-cut decomposition of G.We define contracting a vertex v as identifying 
all vertices of N(v) to a single vertex and deleting v. 

In [25] , Algorithm I A. II is proved to compute a branch decomposition of planar graph G of 
width 21 + 1 in time 0(£n). 

Claim A.l Alaorithm \A.l\ computes a sc- decomposition ofG. 

Proof. Observe that the edge set T* of the dual graph of Rg (the so-called medial graph) 
minus the dual edges of T$ in Line [3] forms a tree due to the acyclicity of T5. Every node of 
T* corresponds to an edge of G and in fact, spans the edge set of G. For turning T* into a 
branch-decomposition we a) bijectively map the leafs of T* to the edges of G and b) make T* 
ternary. For a) we generate for every node v in T* from Line [3H3] one ternary tree (or single 
edge tree) , a local tree C v with v one leaf. In Line [3H3] those local trees are merged from the 
leaves of T* to its root in post order such that each C v contributes to leaf v in the overall 
ternary tree C r . We show now that the such formed branch decomposition actually obeys our 
definition of an sc-decomposition, that is the vertices in each middle set form a cycle in the 
radial graph. Note that every edge of G forms a 4-cycle in Rq- Let e* G E(T*) be the dual edge 

2 In fact the authors construct a carving decomposition out of the spanning tree of the dual graph of the radial 
graph that one obtains after deleting the dual edges of the BFS spanning tree. 
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Algorithm A.l: Computing SC-decomposition. 



Input : Plane graph G, face / G F(G) , radial graph Rq ■ 
Output: Branch-decomposition of G of width at most 21 + 1. 

1 Construct embedded BFS tree T$ of Rq at root /. 

2 Set T* = R* G \ E(T S )* the dual graph of R G without the edge set dual to T s 

3 for every node v in T* do 

4 
5 
6 
7 

8 



if deg(f )t = 1 then 

create C v a single edge, with nodes labeled {v} and {N(v)}; 
else 

create embedded ternary tree C v with |iV(u)| + 1 leaves; 

label one leaf with {v } and the other leaves with N(v) keeping a clockwise order. 
9 (in post order) for every edge {u,v} in T* , where v is the parent node do 
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combine C u and C v by identifying leaf {v } in C u with leaf {u} in C v , and 
contract the identified node and set new tree to be named C v . 



12 Return (C r (for r root of T*)). 



of edge e = {/, g] £ E(Rg \ Tg). Then the union of e and the path through Ts from / over the 
lowest common ancestor of /, g in Tg to g forms a cycle in Rq that separates the two subtrees 
of T* that are separated by e* . Thus, T* already possesses middle sets that form nooses in G. 
However T* is not ternary since it may have maximum node degree 4. The leaves of each local 
tree C v in Line [3] are embedded in the same order as the inverse of their labels, the neighboring 
nodes of v in T* , and thus we keep the same ordering in the overall ternary tree C r . Every 
edge e of C r comes from an edge of one of the local trees C v , and e separates the neighbors 
Ni(v) from -/V^f) where the disjoint union Ni(y),N2(v) form the neighborhood N(v) of v in 
T* . Like this, T* falls apart into two subtrees each bounded by a cycle in Rq formed similarly 
as above by the union of Ni(v) bounding minimal path in Ts and the path through the edges 
of T* induced by iVj [v] , for i = 1 and i = 2 respectively. 
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